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Abstract. The relation between redshift and the CMB temperature, Tqmb{z) = 
To(l + z) is a key prediction of standard cosmology, but is violated in many non- 
standard models. Constraining possible deviations to this law is an effective way to 
test the ACDM paradigm and search for hints of new physics. We present state-of- 
the-art constraints, using both direct and indirect measurements. In particular, we 
point out that in models where photons can be created or destroyed, not only does 
the temperature-redshift relation change, but so does the distance duality relation, 
and these departures from the standard behavioiir are related, providing us with an 
opportunity to improve constraints. We show that current datasets limit possible 
deviations of the form Tcmb{z) = To(l-|-2;)^-^ to be /3 = 0.004±0.016 up to a redshift 
z 3. We also discuss how, with the next generation of space and ground-based 
experiments, these constraints can be improved by more than one order of magnitude. 



1. Introduction 

Cosmology and particle physics are presently experiencing a truly exciting period. On 
the one hand, both have remarkably successful standard models, which are in agreement 
with a plethora of experimental and observational data. On the other hand, there are 
also strong hints that neither of these models is complete, and that new physics may be 
there, within the reach of the next generation of probes. 

There arc three compelling and firmly established observational facts that the 
standard model of particle physics fails to account for: neutrino masses, the existence of 
dark matter, and the size of the baryon asymmetry of the Universe. For each of these. 
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the model makes very specific statements, failing however to reproduce the experimental 
evidence. It is precisely our confidence in the model and our ability to calculate its 
consequences that lead us to the conclusion that it is incomplete, and new phenomena 
must be anticipated. This is, of course, the reason for the LHC project. 

Similarly, the last decade saw the emergence of the so-called concordance model of 
cosmology. This can reproduce all the available observations with only a small number 
of parameters, but also requires that about 96% of the content of the universe is in 
a form that has never been seen in the laboratory (and is only known indirectly from 
its gravitational properties). It is thought that dark matter is a subdominant part of 
this, while the dominant one is an even more mysterious component usually called dark 
energy. 

In this context, it is important to identify laboratory or astrophysical probes 
that can give us more information about the nature and properties of this still 
unknown physics. In this work we will discuss one such probe — the temperature-redshift 
relation — , and lay the foundations for exploring its cosmological implications. 

One of the most precise measurements in cosmology is the intensity spectrum of 
the cosmic microwave background radiation: the COBE-FIRAS experiment revealed 
a very precise black-body spectrum [1]. However, this measurement tells us nothing 
about the behaviour of the Cosmic Microwave Background (CMB) at non-zero redshift. 
If the expansion of the Universe is adiabatic and the CMB spectrum was a black-body 
at the time it originated, this shape will be preserved with its temperature evolving as 
T{z) = To{l + z). This is a robust prediction of standard cosmology, but it is violated in 
many non-standard models, including string theory motivated scenarios where photons 
mix with other particles such as axions (see [2] for a recent review), and those where 
dimensionless couplings like the fine-structure constant vary [3j. 

A few measurements of T{z) already exist, but the currently large uncertainties 
do not allow for strong constraints on the underlying models to be set. However, with 
future datasets this will become a competitive probe. It is therefore timely to discuss 
what these measurements can tell us about the underlying cosmological paradigms. 

At low redshifts, say z < 1, the T{z) relation can be measured via the Sunyaev- 
Zeldovich (SZ) effect towards galaxy clusters. This method was applied to ground-based 
CMB observations [H |5], which demonstrated its potential. With a new generation 
of ground experiments becoming operational and a forthcoming all-sky survey of SZ 
clusters to be carried out by Planck [6j, the potential of this method will come to 
fruition. At higher redshifts, z > 1, T{z) can be evaluated from the analysis of quasar 
absorption line spectra which show atomic and/or ionic fine structure levels excited by 
the photon absorption of the CMB radiation (The CMB is an important source 
of excitation for species with transitions in the sub-millimeter range.) Although the 
suggestion is more than four decades old, measurements (as opposed to upper bounds) 
were only obtained in the last decade, and the best ones so far still have errors at the 
ten percent level 

Here we will study these issues in detail, but we will also place them in a wider 
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context. For example, in models where photons can be created or destroyed, not 
only does the temperature-redshift relation vary, but so does the distance duality 
relation (also known as the Etherington relation [9]), and these two different departures 
from the standard behaviour are quantitatively related. One issue that has been 
overlooked so far is that in such models, where photon number is not conserved, 
this relation between T{z) and distance duality provides us with an opportunity to 
improve constraints. By combining data from different observations one not only reduces 
the statistical uncertainties on underlying phenomeno logical parameters but, given the 
different nature of both observational datasets, one also has a much better control over 
possible systematics. 

We therefore discuss in detail the origin of the above relation, as it can be a unique 
consistency test for the standard paradigm and, at the same time, a valuable tool for 
probing new physics beyond the standard model. We also study further imprints of 
these models in the CMB, and present forecasts for improvements that Planck, as well 
as planned Baryon Acoustic Oscillations (BAO) missions and spectrographs planned 
for the VLT and the E-ELT, will soon make possible. Last but not least, we derive 
the strongest constraints to date on deviations of these relations from their standard 
behaviour, and quantify the improvements to be expected from the aforementioned 
forthcoming experiments. 

2. Theoretical motivation 

There are several examples of non-standard, but theoretically well-motivated, physical 
processes that could affect the cosmological temperature-redshift relation. Constraining 
deviations from the standard law T{z) = Tq{1 + z) therefore provides a invaluable tool 
for probing physical theories. Examples of scenarios that could be constrained include 
decaying vacuum cosmologies/photon injection mechanisms, couplings between photons 
and axion-like particles, modified gravity scenarios, and so on. In this section, we discuss 
two examples of how these models affect the T{z) relation. 

2.1. Adiabatic Photon Injection 

Naively, the most obvious way to violate the T{z) = Tq{1 + z) relation is when there is 
energy injection into the CMB (say from a decaying scalar field) or, conversely, when 
photons are destroyed. The Planckian form of the spectrum is preserved if photon 
creation is adiabatic, that is, if the entropy per photon remains constant. Lima et al. 
Pm [TT] have studied this in the context of decaying vacuum cosmologies; here we review 
some of their formalism. 

Generically, such a process of particle creation can be described by saying that the 
energy density evolves as 

p + 3H{p + p) = C (1) 
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while the particle number density obeys 

h + 3Hn = ^ , (2) 

where C, \l/ are real functions of time. Deviations from particle number conservation 
can then be quantified by the phenomenological parameter 

0<P.^<1, (3) 

which can be time-dependent. 

In general the temperature evolves as 

T f dp\ h ^ 



T \dpj^n nT{dp/dT)^ 



nC 



(4) 



but note that the second term is zero if C and are related by 

C = P±^^, (5) 

n 

or, in terms of the parameter defined in (|3]), C = 3H(p + This corresponds to 
adiabaticity: the specific entropy per particle of the created particles remains constant, 
so new particles are actually created in equilibrium with already existing ones. In this 
case the temperature evolution equation reduces to 

T f dp\ h 



T I ^ (6) 

as in the standard cosmological model. Recall that, in the standard model, both the 
equilibrium relations and the Planckian form of the spectrum for photons are preserved 
in the course of expansion. The former can be seen directly from equation ([6]), while 
the latter is the result of the kinematical condition for FRW geometry, v oc a~^, 
and of photon number conservation, which together imply T oc (away from mass 
thresholds). With adiabatic photon creation, equation ^ remains valid - as we just 
saw - and thus the equilibrium relations are preserved. 

Indeed, if one considers the generic equation of state p = (7 — l)p one can integrate 
the above equations to get 

T oc rC'^ , p oc . (7) 

However, since the entropy per photon is conserved, the temperature now obeys 
Ta 

^ = const, (8) 

where N{z) is the comoving photon number, which changes in time as new photons are 
being injected in equilibrium. It then follows that the dimensionless frequency relevant 
for the SZ effect (to be considered in Q evolves as 



x = Xo{l + z)-^, (9) 
1 r 



-CMB 

with 

1/3 



TcMB=To(l + z) ^ . (10) 
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XmT = 0.289 



Thus, it is now hu{N/NoY^^ /kT in the exponential of the photon distribution function 
that stays constant during cosmological evolution, and so a generalised Planck-type 
spectrum is preserved, which is why the equilibrium relations are still recovered, 
e.g. p (x T'^ for a radiation fluid (where the proportionality factor is the radiation 
density constant). 

Although one cannot distinguish at present the usual Planckian spectrum from 
such a generalised one, this is not necessarily true at higher redshifts. For example, the 
wavelength Am of the peak of the distribution is now 

■^1'". (11) 

which is a generalisation of Wien's law (and naturally reduces to it in the standard 
case). This could in principle be observationally tested. 

We can illustrate these points with two simple examples. Assuming a radiation 
fluid, p = p/3 we have 

f ^ C 

- = -H + — = -H + — (12) 

and for a constant /3 we find 

p + 4(1 - /3)i7p = (13) 

n + 3(1 - /3)iJn = (14) 

T{z) = To{l + zf-^- (15) 

the dimensionless frequency is therefore 

- = {l + zf. (16) 

For a more general equation of state p = (7 — l)p we have 

T oc a-3(^-i)(i-« . (17) 

Thus, for any given redshift the temperature of the expanding universe is slightly lower 
than in the standard case. This is the phenomenological parameterisation that has been 
used by almost all previous authors. Note that for many realistic models this is not an 
accurate description (as (3 is in general time-dependent), but it is adequate at sufficiently 
low redshifts. Since, currently, cluster data probe down to redshifts less than unity with 
sensitivities in temperature of order a few percent or worse (see we will adopt this 
parameterisation for the purposes of the present paper (but will revisit the issue in a 
follow-up publication). As cosmological data improve and extend to higher redshifts 
(see time variation of the parameter (3 will also be constrained. 
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2.2. Photon dimming/absorption & Axion-Photon Couplings 

The effect of photon dimming/absorption can arise in a wide range of astrophysical and 
high-energy physics scenarios, ranging from photon absorption by grey dust to photon 
conversion into a different particle species, like for example photon-to-axion conversion 
in the presence of a magnetic field. 

Grey dust has been invoked as an alternative explanation of the observed dimming 
of Type la Supernovae [12], but it is now understood that it cannot fit high redshift 
data (like for example the Union sample [13]) by itself, i.e. without some contribution 
from a cosmological constant-like fluid [13]. Axion-Like-Particles (ALPs), on the other 
hand, arise in a wide range of well-motivated high-energy physics scenarios, including 
string theoretic models where ALPs appear as zero modes of various antisymmetric form 
fields [in]. Like dust, they can lead to dimming of Type la Supernovae [TB], but this 
cosmological scenario is also strongly constrained [T71 [T8t [T9] . For a recent review of 
ALPs and related laboratory, astrophysical and cosmological constraints, see [2]. 

Irrespectively of its microphysical origin, photon dimming violates photon number 
conservation, and so it can be described macroscopically in a very similar way to photon 
injection, by simply allowing the parameter [5 in equation ([s]) to be negative. In 
particular, the balance ( [l4| ) and temperature evolution ^ equations are valid with 
constant negative /3 (note that, as before, constant (3 is only an approximation) and the 
temperature-redshift relation is again: 



Crucially for this work, violations of photon number conservation also give rise 
to deviations from the standard relation [9] between luminosity and angular diameter 
distance, which can be independently constrained [201 [211 1221 IH HH]- With current 
data, it is sufficient to parameterise the luminosity- angular diameter distance relation 
as: 



where e is constant. The value e = corresponds to the standard di — dA law. Since e 
and (3 are related by the underlying physical model but can be constrained independently 
through measurements of very different systematics, this provides a very promising tool 
for carrying out a consistency test of the standard cosmological scenario and constraining 
physics beyond the standard model. However, such complementary measurements often 
cover a very different wavelength range so one must be careful when comparing direct 
and indirect constraints. For example, if one uses Supernova (SN) data to constrain the 
parameter e at optical wavelengths within a model of chromatic axion-photon mixing, 
one should take into account the wavelength dependence of the model when translating 
this bound into a constraint in the CMB temperature-redshift relation. In general, one 
may expect photon-dimming to be stronger at high photon energies, so indirect bounds 
on T{z) coming from SN brightness measurements can be assumed to be conservative. 
In [|5] we will present the constraints obtained from such an analysis with current data. 






(19) 
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and in ^ we examine the prospects of improving these constraints by future SN and 
BAO data. 

In passing, we note that another way to see that the temperature-redshift relation 
becomes modified in models violating photon number conservation is by considering the 
more obvious effect of the distance-duality violation in these models. The source of this 
distance-duality violation is the difference between the true and observed (dimmed) 
photon flux, which tricks one to infer a larger luminosity distance. As flux scales 
inversely with the luminosity distance squared, this change in the inferred distance 
would also affect the CMB fiux-redshift relation, and so the CMB temperature as a 
function of redshift. 

Before moving to the study of the various constraints imposed on (3 and e, let us 
explore their relation in more detail. Consider the CMB spectrum in the presence of a 
dimming agent, like for example photon absorption due to a dust field or conversion into 
ALPs. Imagine that the spectrum is Planckian at the epoch T so that the the number 
of photons per unit volume per frequency interval is: 

Stt u'^du 

^(^,)rf^,= (20) 

If photon number was conserved, then at a later epoch T' the number density per 
frequency interval, n'{v')dv\ would be related to (20) through the volume rescaling 
{a/a'Y; this can be absorbed as frequency dependence, z/' = [ajd^v, yielding the 
standard result of a Planckian spectrum niv'^dv' of temperature T' = {a/a')T. With 
dimming, however, photons can be lost in flight and the spectrum will generally be 
distorted. Introducing a photon survival fraction /surv ^ 1 between the epochs T and 
T' we now have: 

n\p')dp' = f surv n{iy)diy {a/a'f = ^^j^rj^fr^ < n{p')dp' , 

where, as before, z/' = v{a/a') and T' = {a/a')T. Thus the spectrum is distorted. 
However, if the photon survival fraction does not depend on frequency (of course, 
it can — and generally will — depend on redshift), then one can change variables to 
i^" = ifsurvY^^ i^' , yielding: 

n'{iy')diy' = n{iy")du" , 

that is, a Planckian spectrum of temperature T" = {fsurvY^^T' = {a/ a'){fsnrvY^^ T . 
This is the 'generalised' Planck distribution of reference jTU], discussed in the previous 
section. In fact, remembering that T' corresponds to an epoch later than T, we have: 

T{z) = {U^„Y^/\l + z)T^. (21) 

The relation between dimming and temperature distortion then becomes clear. In 
the parameterisation (15), we have {fsunY^^ = (1 + ■, with (3 negative. On the 
other hand, the photon survival probability /surv enters linearly the luminosity, so the 
luminosity distance scales as (/surv)~^^^- Therefore, from equation (19) it follows that the 
photon survival probability corresponds, in the e-parameterisation, to /surv = (1 + ^)"^". 



Constraints on the CMB temperature-redshift dependence 



8 



Thus, the relation between the parameter /3 of equation (15) and e of references [Ml [T9] 
is simply: 

/3 = -^e. (22) 

Note that in the above discussion we have assumed that the photon survival 
probability is independent of wavelength, which guaranteed a simple thermal spectrum. 
Deviations of the CMB from a thermal spectrum have been constrained down to the 10~^ 
level by the COBE-FIRAS measurements, both in the case of a Comptonised spectrum 
and of a Rose- Einstein spectrum with non-trivial chemical potential [H |23] • This dataset 
has also been used recently to constrain frequency-dependent dimming in the context 
of photon- axion couplings p!H]. Finally, we have assumed a simple redshift dependence, 
parameterised by equation (19) with constant e. Again, like in the case of constant (3 
(section [2T ), this is justified given current error bars in distance determination and the 
limited redshift range currently covered, < z < 2 [13]. Forecasts for future data will 
be discussed in ^ 



3. Constraints from quasar absorption line spectra 

The most accurate value of the local CMB temperature measured by the CORE 
experiment is 

TcMB = 2.725 ±0.002 is:, z = 0. (23) 

Additional local but extrasolar values of the background radiation temperature can be 
estimated from observations of interstellar molecular clouds [23] and of the Magellanic 
Clouds [25j. These extrasolar and extragalactic values are in good agreement with 
the COBE estimation. Cosmological models predicting non-Planckian spectra of the 
background radiation at z = can be ruled out if they deviate by more than about 
1% from the blackbody spectrum. However, these local observations do not allow one 
to distinguish between the standard model and cosmological models with a blackbody 
spectrum but different T{z) dependences. 

At low redshifts (0 < z < 1), the functional scaling of the CMB temperature can be 
estimated from measurements of the Sunyaev-Zeldovich effect in the direction to clusters 
of galaxies, as will be discussed in the next section. At higher redshifts {z > 1), the 
CMB temperature can be evaluated from the analysis of quasar absorption line spectra 
which show atomic and/or ionic fine structure levels excited by the photo-absorption 
of the CMB radiation. This is an important source of excitation for those species with 
transitions in the sub-millimetre range. This is the case for atomic species whose ground 
state splits into several fine-structure levels (and also for molecules that can be excited 
in their rotational levels, see below). If the relative level populations are thermalised 
by the CMB, then the excitation temperature gives the temperature of the black-body 
radiation. 

It has long been proposed to measure the relative populations of such atomic levels 
in quasar absorption lines to derive Tcmb at high redshift [2S] • The most suitable are the 
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fine structure levels of the ground states of CI (3P0, 1, 2) and CII (2P1/2, 3, 2) showing 
an energy separation between 24 K and 91 K. The CI lines were used to obtain several 
upper limits, until advances in instrumentation and analysis techniques allowed actual 
measurements, starting with Srianand et al. [7] who obtained 

Te = 10 ±4:K, z = 2.338. (24) 

To distinguish the contribution to the relative population of the fine-structure levels of 
the ground state of CI from competing excitation processes the independent analysis of 
the molecular hydrogen UV absorption lines is important; H2 transitions from different 
low rotational levels may be used to infer the UV radiation field and the gas density in 
the CI-H2 absorbers. These techniques have so far allowed measurements to be made 
beyond z = 3; the earliest currently available measurement is that of Molaro et al. [27] 

T^ = 12.ltl-lK , 2 = 3.025. (25) 

Molecular rotational transitions may also be used to constrain the cosmological 
temperature-redshift law. So far, absorption and emission lines of CO, OH, CS, HCN, 
HCO^ , H2O, NH2, and other molecules have been observed in distant galaxies and 
quasars up to z = 6.42 (see [28] for a review). The absorption line observations are most 
suitable to the radiation temperature estimates since, as a rule, molecular absorption 
arises in the gas components with the lowest kinetic temperatures. Srianand et al. [21] 
reported the first detection of CO in a high-redshift damped Lyman-a system, while also 
detecting H2 and HD molecules. The CO rotational excitation temperatures are higher 
than those measured in our Galactic ISM for similar kinetic temperature and density. 
Using the CI fine-structure absorption lines, they show that this is a consequence of the 
excitation being dominated by radiative pumping by the cosmic microwave background 
radiation, and from the CO excitation temperatures they derive 

Te = 9.15 ± 0.72 K , ;z = 2.418. (26) 

Finally, Noterdaeme et al. [8\ have recently reported on a sample of five CO 
absorption systems where the CMB temperature has been measured. We refer the 
reader to this work for further details as well as for some more on the history of these 
measurements and an up-to-date list of all the available ones. They also used their 
sample, in combination with measurements from the SZ effect, to place constraints on 
the phenomenological parameter /3 — we will describe this in the next section. 

4. Constraints from the SZ effect towards clusters 

The possibility of determining TcMsiz) from measurements of the Sunyaev-Zeldovich 
effect has been suggested long ago [201 EI] • The effect - Compton scattering of the CMB 
by hot intracluster (IC) gas - is a small change of the CMB spectral intensity, Al, which 
depends on the integrated IC gas pressure along the line of sight to the cluster. The 
steep frequency dependence of the change in the CMB spectral intensity, Al, due to the 
SZ effect allows the CMB temperature to be estimated at the redshift of the cluster. 
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The differential SZ signal may be written at the redshift of the cluster, including 
relativistic corrections as: 

A/(z) = ^^^^^^^'^^\f'\y r mx) - vjc + Rix, 9, vjc)] , (27) 

where AI{z) is the brightness change between the centre of the cluster and blank sky, 
as measured at redshift z, t is the optical depth, Tcmb(-2) is the CMB temperature at 
redshift z, x = kBTcMsiz) ' ^ ~ kBTe/nieC^ with electron cluster temperature, Vz the 
radial component of the peculiar velocity of the cluster, and the R{x,9,Vz/c) function 
includes relativistic corrections. 

If we assume that Tcmb scales with z as Tcmb{z) = Tcmb(0)(1 + zY''^, while the 
frequency scales as (1 + z) as usual, then A/(z) scales like (1 + z)^^'^"^'^ and we obtain, 
at the level of the solar system 

^^(0) ^ 2ikTc^P^_^^^^ _ ^^^^ ^ ^^^^^^ ^28) 

where x' = and T^j^^ = Tcmb(0)(1 + -z)^^ will be slightly different from the 

local temperature Tcmb(O) as measured by COBE. In this way it is possible to measure 
the temperature of the CMB at the redshift of the cluster, thus directly constraining 
scenarios like those discussed in the previous section. 

Let us consider the bound on /3 coming from Tcmb(^) constraints, at redshifts in the 
range z = 0.023 — 0.546, from multi-frequency measurements of the SZ effect towards 
the 13 clusters of Ref |5]. By fitting the Tcmb (-2^) data points and relaxing the condition 
of a positive prior on /3, that is allowing for photon dimming/absorption (corresponding 
to /3 < 0) as well as photon creation (/3 > 0), we get /3 = 0.065 ± 0.080. 

By adding these measurements to higher-redshift ones coming from spectroscopic 
measurements (discussed in the previous section) involving atomic carbon and CO 
absorption lines along the line of sight of quasars, Ref. [H] subsequently improved the 
constraints on /3, obtaining 

/3 = -0.007 ±0.027. (29) 

Despite the higher precision of Tcmb (2;) measurements from the SZ method with 
respect to that of the quasar absorption lines method, the present improvement on 
the constraints on the parameter /3 depends mainly on the higher lever arm due to the 
exploration of the distant universe, thanks to the observation of high redshift absorbers. 



As we will show in ^6.1 forthcoming datasets will soon lead to tighter constraints. 



5. Constraints from distance measurements 

As briefly alluded to in [|2| an indirect way for constraining the temperature-redshift 
dependence is through the study of possible deviations from the duality relation [9] 

I Repeating their fit to parameter /3 we find an uncertainty of 0.028 ratlrer than 0.027; this is likely due 
to round-off errors in the temperature measurements used as input. The difference is insignificant for 
the purposes of our subsequent analysis, but for the sake of consistency with the pipeline for comparing 



to future measurements we will use 0.028 in Figs. 10 and 11 
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between luminosity and angular diameter distance: 

dLiz) = il + zYdAiz). (30) 

This equation holds for general metric theories of gravity, where photons travel along 
unique null geodesies, as long as photon number conservation and local Lorentz 
invariance are respected. Since Lorentz violations are strongly constrained at low 
energies [32], and in particular at optical wavelengths, the determination of di from SN 
observations can be used |20l [211 l22l [lH HI] together with other distance measurements 
in order to place direct bounds on photon number violation through Eq. (30 ). Physically, 
such a violation can arise from photon absorption, e.g. grey dust [12j, or from more 
exotic effects like photon conversion into axions \16\ . 

Systematic violations of (30) give rise to an apparent opacity effect in the observed 
luminosity distance. Indeed, if photons were lost along the line of sight, the inferred 
luminosity distance would be related to the true one through a multiplicative factor: 

Note that the 'opacity' t{z) can be negative, allowing for apparent brightening of the 
source, as would be the case, for example, if exotic particles were also emitted from the 
source and later converted into photons [33] . 

As discussed in ^ photon number violation would also give rise to a corresponding 
distortion of the photon temperature-redshift relation, and this allows us to combine 
different observational probes to constrain such models. 

In references [TH [19] the authors used Type la SN data (specifically the Union 
2008 dataset [13]) in combination with measurements of cosmic expansion H{z) from 
differential ageing of luminous red galaxies [311 EHl [36], and obtained constraints on 
opacity up to 2; ~ 2 through equations (30 31). In the simplest case, one can adopt the 
parameterisation: 

dL{z) = dA{z){l + zf+% (32) 



allowing violations of Eq. (30) through a single parameter, e. For low redshifts, 
this parameter can be directly translated into an opacity function t{z) as it simply 
corresponds to the first term in the Taylor expansion t{z) = 2ez + Olez"^). Different 
parameterisations were also considered in ^9j, corresponding to specific theoretical 
models of opacity such as photons mixing with massless axion-like particles, chameleons 
and mini-charged particles. 

If an opacity source like photon-axion mixing affects SN observations, it should 
also have an impact on CMB photons. From the above discussion (and that of ^2.2) 
it is then clear that constraints like those in [TH [T^] can be used to place indirect 
bounds on possible deviations from the standard temperature-redshift law, as expressed, 
for example, in Eq. (15). In particular, if opacity is achromatic, which can be the 
case, for example, for photon-axion mixing in sufficiently low intergalactic electron 
densities [IQ], the parameterisation (32) corresponds simply to (3 = — 2e/3, as we saw 
in ^2.2 More generally, however, the effective opacity can be expected to be smaller 
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at low energies, so that distance duality violation bounds yield conservative constraints 
for the allowed CMB temperature-redshift distortion. There are also other bounds 
on photon-axion mixing, coming from constraining spectral distortions of the CMB 
radiation yL8j. Note that for CMB photons to be converted into axions, a background 
magnetic field is required at the corresponding redshift. Therefore, the bound on the 
parameter e constrains a combination of the background magnetic field and the axion- 
photon coupling [T71 UHl US] (and at the same time, through its relation (22) to /3 it 
places an indirect constraint on the variation of T{z) as discussed above). Intergalactic 
magnetic fields have not been directly observed but there are hints that they exist [37] . 
For further discussion and relevant constraints, see [3Sl ES]- Finally, note that axions 
can also convert back into photons - a second order effect which could be significant for 
large mixing probability. 

Let us consider the constraint on /3 coming from the distance duality bounds of 
Refs. [lll[T9]. With (3 = — 2e/3, combining the most recent SN and H{z) data (namely, 
the SCP Union2 Compilation [H] together with the latest H{z) [36] data and Hubble 
parameter determination [12]) yields 

/3 = 0.01 ±0.04, (33) 

at 95% confidence. 

This has been obtained by considering fiat ACDM models and marginalising over 
the matter density Qrn and Hubble parameter Hq. Fig. [T] shows the relevant constraint 
on the P — Qm plane after marginalising over Hq (left) and the constraint (33) on /3, 
obtained by marginalising over both Hq and Q^- On the left, the dark blue contours 
are the 1 and 2-a joint (2-parameter) confidence levels for the SN data, the lighter blue 
regions show the corresponding constraints from H{z) data, and the solid black lines 
show the combined SN+if(2;) constraints. Again, as we will discuss in ^6.2, there are 
very good prospects for future improvements. 



6. Forecasts 



Having discussed the current constraints on the parameter /3, we now move on to study 
in detail the prospects for improvements coming from the next generation of space and 
ground-based experiments. We will discuss in succession three different probes, which 
together span the redshift range z ~ 0—4. By combining data from different observations 
one can therefore explore thoroughly a wide range of redshifts and significantly reduce 
the statistical uncertainties on the underlying phenomenological parameters. Moreover, 
given the very different nature of each type of measurements, one also has a much 
better control over possible systematics — this is particularly important in the event of 
a detection of deviations from the canonical behaviour. 
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Figure 1. Constraints from SN-|-i?(z) on the parameter f3, parameterising violations 
of the temperature-redshift relation as T{z) = To(l -I- z)^'^. Left: Two-parameter 
constraints on the /3 — plane. Dark blue contours correspond to 68% and 95% 
confidence levels obtained from SN data alone, light blue contours are for H{z) data, 
and solid line transparent contours show the joint SN+H{z) constraint. Right: One- 
parameter joint constraints on /3 marginalised over cosmological parameters. The solid 
line is for SN+H{z) while the dotted one for SN data only. The dashed line shows the 
Ax^ = 4 level. 

6.1. Low redshifts: Planck HFI (clusters) 

Planck HFI |13] was specifically designed from the beginning to measure the SZ effect 
in galaxy clusters [H]: the spectral coverage allows one to explore the positive and 
negative part of the spectral distortion, and it is optimally suited for cluster detection 
and to break cluster parameters degeneracy. 

The full sky survey will provide us with a SZ catalog of thousands of clusters. The 
Planck early results already provide us with a sample of 189 cluster candidates. 
Here we focus on a survey dedicated to a sample of well known clusters, for which X-ray 
and optical information is available, so dealing with a subsample of the final Planck 
cluster catalog (ESZ). 

A catalog of 166 clusters has been built by using BAX (X-ray Clusters Database) 

m- 
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Cluster name 
RA (J2000) 
DEC (J2000) 
z 

Fx 

Reference-Fx 
Lx 

Reference-Lx 
Band-Inf (keV) 
Band-Sup (keV) 
Tx 

O-Tx 

Reference-Tx 
Instrument 

R-core 

Reference-fZcore 



Right ascension 

Declination 

redshift 

unabsorbed X-ray flux in ROSAT band (0.1- 
2.4) keV in units of (10~^^erg/s/cm^) 

X-ray luminosity in the ROSAT band (0.1- 
2.4) keV in 10^^ ergs s'^ 



X-ray gas temperature in keV 



Core Radius (arcsec) 



slope of the gas density profile derived from 
the /3-model fitting 



To avoid confusion, note that this /3-model does not refer to the same (3 parameter 
as is discussed in the rest of the paper. 

For each cluster we derive the following parameters: 



Ueo central electronic density, assuming a 

isothermal /3-model and following [IB] 
Uth central Comptonisation parameter 

Tth central optical depth 

Yint Comptonisation parameter integrated over 

the cluster extent 
Daszx Angular distance 



We have simulated the observations of a sample of about 40 well known clusters, 
already observed by Planck [H], taking into account the Planck HFI instrumental 
characteristics and observing strategy [47]. The frequency bands considered in the 
simulation are the four at lower frequencies (100, 143, 217, 353 GHz), since the remaining 
two higher frequency bands (545, 857 GHz) are best suited for foregrounds extraction 
and are not useful for the reconstruction of the SZ signal. CMB and foregrounds are 
assumed as previously removed. 
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The forecast for the SZ effect signal for the clusters has been obtained from the 
measured X-ray properties, assuming an isothermal model. The ICM pressure profile has 
historically been described by an isothermal /3-model |19]. Recent X-ray observations 
have shown that a /3-profile for gas density is a poor approximation for the cluster's 
profile at large radii, leading several authors to propose more realistic profiles [HUl EI]- 
The use of the /9-model in this work is for consistency with the X-ray derived parameters 
collected in the catalogue, which almost invariably were based on a /3 profile assumption 
for gas. Nevertheless, the analysis procedure to determine T(z) makes use of the central 
T values, and it will not be affected by a different adopted profile. 

We have used the Planck noise model (NET values as reported in |17j) to estimate 
the errors in the observed spectra. The integration time for the mock observation is 
obtained assuming uniform sky coverage and two years of observation. For each cluster 
the observation time is 10s. The dilution effect has been taken into account to estimate 
the error on the SZ signal for each channel. 

The mock dataset was then analysed to recover the original input parameters of 
the cluster. The analysis has been performed through a Monte Carlo Markov Chain 
(MCMC) algorithm which allows us to explore the full space of the cluster parameters 
(optical depth r, peculiar velocity Vpec, electron temperature Te) and the CMB brightness 
temperature at the redshift of the cluster. In the analysis we allowed for calibration 
uncertainty, considered as an uncertain scale factor, which was modelled as a Gaussian 
with mean 1 and 0.1% standard deviatior][§] 

The MCMC generates random sequences of parameters, which simulate posterior 
distributions for all parameters [51] . The sampling approach we used is the one proposed 
by Metropolis and Hastings . Convergence and mixing of the MCMC runs was 

tested through the Gelman- Rubin test [57]. We included a prior over the cluster gas 
temperature, as provided by X-ray data. For the radial component of the cluster peculiar 
velocities the prior is a theoretical one, a Gaussian with a universal vanishing mean and 
with a 1000 km/s standard deviation. Clusters with almost flat r posterior are excluded 
from the sample (this is the reason for which the sample is different when the kinematic 
component is included from that when it is not included). 

As was already noted in [3] there is a degeneracy between Tcmb(-2) and Vpec, see Fig. 
[2| In order to reduce the impact of this degeneracy and then to reduce the uncertainty 
in the determination of Tcmb(-2^), a better knowledge of the peculiar velocity is required 
or it is necessary to remove the kinematic component from the thermal component, 
together with the CMB intrinsic anisotropy. With only Planck measurements, given the 
multifrequency coverage, we can envisage the second option. Kinematic SZ (KSZ) and 
the intrinsic CMB anisotropy have the same spectral shape in the non relativistic limit; 
to separate the two components in nearby clusters or in distant clusters using beams 

§ The adopted value for the absolute calibration accuracy is in a way too optimistic with respect to 
HFI early maps (< 2%) [S^, but still very conservative with respect to expected performances |53j . 
However, the impact of absolute calibration accuracy will mainly influence the r parameter estimation, 
since A/ depends linearly on r at the first order. 
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Figure 2. Parameter correlations for a simulated cluster: the contours show the 
68% and 95% confidence limits fi-om the marginalised distributions. Black contours 
are obtained allowing for a peculiar velocity prior with vanishing mean and standard 
deviation 1000 km/s; cyan contours are obtained allowing for a peculiar velocity 
prior with vanishing mean and standard deviation 100 km/s. Tcmb corresponds to 

r(z)/(i + z) 



larger than few arcmin it is necessary to rely on the very small spectral distortions of 
the kinematic SZ due to relativistic effects [58]. For Planck it is safe to consider that 
in a large majority of cases the spectral distortions of KSZ due to relativistic effects are 
too small with respect to CMB confusion and noise level to allow disentangling KSZ 
and CMB intrinsic anisotropy, i.e. the kinematic component is removed as a first step 
in cleaning maps of CMB contribution. 

Nevertheless, in the following we will consider both cases, thus fitting either four 
parameters (r, Tg, Vpec, Tcmb(^)) or three parameters (r, Te, Tcmb{z)). 

In Figs. [3] and |4] we present the results of the parameter estimation analysis for a 
single cluster, both with and without the kinematic component. In Figs. |5] and [6] we 
show the residuals for cluster parameters for the whole sample. 
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Figure 3. An example of parameter extraction for a single cluster, with kinematic 
component included. 
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Figure 4. An example of parameter extraction for the same cluster of Fig. |3] with 
kinematic component not included. 
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Figure 5. Residuals for 32 clusters, with kinematic component included. All the 
clusters with almost flat r posterior have been excluded from the original 42-cluster 
sample. 
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Figure 6. Residuals for 37 clusters, with kinematic component previously removed. 
All the clusters with almost flat r posterior are excluded from the original 42-cluster 
sample. 

To obtain (3 we have performed a fit of the T(z) data points (see Figs. [7]and|8]). The 
final P value we get by fitting the T{z) data points obtained with the MCMC treatment 
is (3 = —0.047 ± 0.079 when the kinematic component is included (from a final sample 
of 19 clusters, selected with the condition of non fiat r posterior and S/N > 6) and 
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Figure 7. Left: Tcmb vs z, with tlie kinematic component included. Right: Posterior 
of the /3 parameter, as obtained by performing a fit of the T{z) data points. 



5.0 




Figure 8. Left: Tcmb vs z, with the kinematic component not inchided. Right: 
Posterior of the /3 parameter, as obtained by performing a fit of the T{z) data points. 



(3 = —0.003 ± 0.016 when the kinematic component is previously removed (from a final 
sample of 37 clusters, selected with the condition of non fiat r posterior and S/N > 6). 
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In conclusion, if the kinematic component is removed altogether with the CMB 
primary anisotropy component, with Planck we can reach 0.6% sensitivity on Tcmb{z) 
measurements; otherwise the sensitivity will be around 7%. With only tens of clusters 
we can in principle get better constraints on /3 than the current results with SZ+Atomic 
carbon+CO. 

6.2. Intermediate redshifts: EUCLID/SNAP (distance measurements) 

In this section, we show forecast constraints on the temperature-redshift relation (in 
particular on the parameter /3 of which could be achieved by combining H(z) 
measurements from upcoming spectroscopic BAO surveys with future SN data. The 
next decade will see a dramatic improvement on H{z) and angular diameter distance 
data at redshifts z <2, notably through ongoing and upcoming BAO surveys like BOSS 
[nn] and EUCLID [20] ■ Similarly, future SN missions (e.g. SNAP [61J) will dramatically 
reduce the errors in SN brightness data. 

As discussed in [19], BOSS will not significantly improve opacity bounds (on which 
our constraint on (3 are based) with respect to current H{z) cosmic chronometer data 
(i.e. from differential ageing of luminous red galaxies |3^), because it will be restricted 
to redshifts z < 0.7. 

On the other hand EUCLID - a combination of the earlier SPACE [62] and DUNE 
[63] missions - will reach much higher redshifts and is expected to dramatically improve 
these constraints. Aiming for launch in 2019, it would cover about 20,000 deg^ of sky 
providing around 150 million redshifts in the range z < 2. Here, we consider forecast 
constraints from EUCLID and a Supernova SNAP-like survey (or dark energy task force 
stage IV SNe mission) [6^ . 

We use the code developed by Seo & Eisenstein [HS] to estimate the errors in radial 
distances achievable by using BAO as a standard ruler. Fig. [9] shows our forecasted 
constraints on the parameter /3, using modelled BAO data with forecasted errors for 
EUCLID, combined with modelled SN data and errors for a SNAP-like survey. 

On the left panel, light blue contours show the 1-cr and 2-a (2-parameter) 
constraints on the plane from EUCLID only, darker blue contours show the 

corresponding constraints from SNAP, while solid line transparent contours show the 
(2-parameter) joint EUCLID-(-SNAP forecast constraints. To make a more direct 
comparison we have also shown the corresponding constraints obtained from current 
data, namely 'cosmic chronometer' H{z) (dashed), SN (dotted), and joint H{z)+S'N 
(dot-dashed), discussed in ^ The right panel shows the relevant 1-parameter constraint 
on P, after marginalising over Qm, for EUCLID-I-SNAP (solid) and current if(z)+SN 
data (dashed). The dotted line is the 95% confidence level, Ax^ = 4. 

Overall, the improvement with respect to current constraints by combining 
EUCLID+SNAP will be quite significant, with the area of the joint constraint in Fig. |9] 
(left) reduced by a factor of a few decades. The one parameter constraint on /3 will be 
improved by a factor > 5, reaching | A/3| ~ 0.008 at 95% confidence. This is competitive 
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Figure 9. Future constraints on tiie parameter f3 from EUCLID and SNAP. Left: 
Two-parameter constraints on the /3 — plane. Dark blue contours correspond to 
68% and 95% confidence levels from SNAP alone, light blue contours are for EUCLID, 
and solid line transparent contours show the joint SNAP-I-EUCLID forecast constraint. 
Also shown are current constraints from H{z) 'chronometer' data (dashed), SN data 
(dotted), and joint i7(z)+SN (dot-dashed), presented in ^js] Fig. [l] Right: One- 
parameter joint constraints on /? marginalised over cosmological parameters. The 
solid line shows the forecast constraint from EUCLID-I-SNAP, while the dashed line 
corresponds to the current constraint from iJ(2;)-|-SN, discussed in ^ Fig. [I] The 
dashed line is the 95% confidence level. 



to the Planck HFI result, discussed above. 

6.3. High redshifts: ESPRESSO and CODEX (spectroscopy) 

ESPRESSOU (for the VLT) and CODEXU (for the E-ELT) are two forthcoming 
ESO high-resolution, ultra-stable spectrographs. Although their common cosmology- 
related science driver is the precise spectroscopic measurement of nature's fundamental 
couplings (particularly the fine-structure constant a and the proton-to-electron mass 
ratio /i, see |66]), they will be in a unique position to carry out precise measurements 
of T{z) at high redshift. 

As discussed in ^ there are currently 5 CO absorption systems, in the redshift 
range z ~ 1.5 — 3.0, where the CMB temperature can be measured with an uncertainty 

ATNow~0.7ir. (34) 

Based on current plans for ESPRESSO and CODEX [HTl EHl ES], we can estimate that 

II See http://espresso.astro.up.pt/' 

^ See http://www.iac.es/proyecto/codex/ 
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Figure 10. Forecasts for direct constraints on the phenomenological parameter /3 
for ESPRESSO and CODEX, compared witlr the current uncertainty (coming from 
combining all available direct measurements, cf. SQ. 

they will be able to reduce this uncertainty to, respectively, 

ATesp 0.35 K (35) 

and 

ATcoD 0.07 K . (36) 

Given the planned redshift range of both spectrographs, their measurements will 
on average be done at higher redshifts, and here we will assume a redshift range 
z ~ 2.8 — 4.0; from a theoretical point of view, going to higher redshifts is obviously 
desirable since they provide a bigger lever arm; however, one also has to keep in mind 
that these systems will be fainter. 

Having said that, it is important to realise that the bottleneck here is not the 
amount of telescope time required to observe these systems (although that naturally 
grows as systems become fainter). In fact, in some systems that are observed with the 
aim of measuring n one can also measure T{z), so there could in fact be no extra cost in 
terms of telescope time. Instead, the bottleneck is simply finding more systems where 
these measurements can be made. Ongoing surveys such as SDSS-III BOSS [59] can 
play an important role in this endeavour. For the purposes of the present analysis we 
will assume that ESPRESSO will accurately measure T[z) in 10 systems, while CODEX 
will measure 20 systems. 

One can then generate mock catalogs of T{z) measurements, assuming the standard 
scenario as a fiducial model and the above temperature uncertainties and redshift ranges, 
and determine the constraint on /3 that such a catalog can yield. For simplicity we also 
assume a uniform probability in the redshift distribution of the sources. We will also 
assume that the uncertainty in the measurement at z = remains unchanged. Naturally, 
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the constraint will have a mild dependence on where in the allowed redshift range the 
small number of sources happen to fall, but by generating a large number of realisations 
one can infer a representative uncertainty on /3. 



The results of this analysis are summarised in Fig. [TOj ESPRESSO is expected to 
improve on the constraints coming from currently available spectroscopic measurements 
by a factor of 3, and CODEX should improve on ESPRESSO by another factor of 3. 
As stated above, if more systems are found where the spectrographs can make these 
measurements, one could in principle further improve the sensitivity on f3. 

7. Summary: current and future constraints 

We are now in a position to summarise the constraints on allowed deviations from 
the standard temperature-redshift relation. By combining the direct constraints from 
Noterdaeme et al. ^ with the indirect ones that we have obtained in ^ we finally 
obtain the weighted mean result 

/3 = 0.004 ±0.016, (37) 

which is a 40% improvement on the direct constraint. We note that the three 
observational methods (clusters, distance measurements and spectroscopy) are nicely 
complementary, not only in terms of possible systematic uncertainties but also in terms 



of the redshift ranges covered. These results are summarised in Fig. 11, and compared 
with the forecasts for the various future experiments that we discussed in the previous 
section. 

For experiments whose results will be available before 2020, we have considered 
Planck HFI and the ESPRESSO spectrograph. One can see that Planck alone should 
be able to do as well as all the current data, with only tens of clusters. When ESPRESSO 
becomes available, it will allow a gain of a factor of 3 in sensitivity relative to the current 
direct constraints, and a factor of 2 relative to the constraints available just before it. 

We point out that the Early SZ Planck catalogue [TD] consists of 189 clusters and it 
is obtained with a selection criterion of S/N>6, thus implying that the SZ cluster sample 
for which Tcmb can be extracted is at least a factor of 4 larger than the one presented 
in this work; assuming that the spanned redshift range is the same as the one taken 
into account here (thus sensitivity on Tqmb{.z) measurements is almost unchanged) then 
we have a 50% improvement on the constrained /3 (o"/3 ~ 0.008). In principle, spatially 
resolved spectroscopic observations of galaxy clusters (as proposed with SAGACE [71] 
or MillimetroE[^ would allow to further improve these constraints. Finally, we note that 
there are also other techniques that can in principle be used to obtain these constraints 
from SZ clusters [72] . 

In the longer (post-2020) term, CODEX and a combination of EUCLID and (some 
form of) SNAP will bring significant further improvements. Here the order in which they 
will become available is uncertain (as are, to some extent, their detailed characteristics). 

See httpi/www.sron.rug.nl/millimetron 
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Figure 11. Comparing tlie current one-sigma uncertainties on the parameter /3 (from 
direct measurements and from a combination of direct and indirect ones) with that 
achievable by ongoing and future experiments, specificaUy Planck HFI and ESPRESSO 
in the near future and EUCLID/SNAP and CODEX in the longer term. 



but their combined results are expected to bring a gain of about an order of magnitude 
relative to current sensitivities and of factors between 2 and 4 relative to the constraints 
available just before them. 

Our analysis may in some ways be too simplistic, but we emphasise that in 
other ways it is fairly conservative. This is particularly the case when it comes to 
the assumptions on the number of systems in which the measurements can be made. 
The Early SZ Planck catalogue is evidence of this fact for the low redshift range, 
but the same is true for the spectroscopic measurements at high redshift. Given the 
exquisite resolution and stability of the forthcoming ESO spectrographs and the large 
redshift lever arm they can probe, the most cost-effective way to further improve these 
constraints is undoubtedly to identify further systems where ESPRESSO and CODEX 
can make these measurements. 

8. Conclusions 

In this paper we have explored novel techniques for constraining physics beyond the 
standard model, focusing, in particular, on cosmological scenarios that can violate 
photon number conservation. These include - but are not limited to - models in 
which photons mix with axion-like particles, decaying vacuum cosmologies and other 
photon injection mechanisms, models with astrophysical dust, and so on. By noting 
that such models can simultaneously modify the cosmological distance duahty relation 
and the CMB temperature scaling law (these modifications being parametrically related 
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to each other), we have initiated a programme of studying the consistency of these 
models through a combination of direct and indirect probes. Relevant probes include 
SN brightness measurements (yielding luminosity distances), galaxy ageing and/or 
BAO techniques (giving radial and angular diameter distances), SZ measurements of 
galaxy clusters (providing T{z) at low redshifts < z <1) and quasar absorption hne 
spectroscopy (measuring T{z) at higher redshifts up to 2; ~ a few). 

This significantly enlarges ones' toolbox for studying cosmological models beyond 
the standard paradigm and leads to a notable improvement on current constraints 
on such models. Indeed, the probes used are complementary, each having different 
systematics and/or redshift cover, so combining data from several probes allows one to 
reduce systematic uncertainties and obtain more stringent consistency checks, as well as 
improved constraints on the models under study. The combined bound (37), which we 
have obtained on the parameter /3 quantifying deviations from the standard T{z) law 
(see equation (jlsj)), is a 40% improvement over the corresponding direct constraints in 

The potential of this programme is enormous. Considering ongoing and future 
missions (Planck HFI, EUCLID, SNAP, ESPRESSO, CODEX) we have obtained 
forecast constraints/errors on the parameter /3. The expected improvement in current 
1-a error bars can be up to two orders of magnitude, as summarised in Fig. 11 There 
is also room for expanding further the current toolbox by accommodating more probes, 
again with different systematics and redshift cover. For example, the position of the 
first acoustic peak in the CMB is also sensitive to distance duality violations and so 
provides a new tool at the redshift of last scattering [73] . 

In this work, we have only tried to demonstrate the effectiveness of these techniques 
by considering simple parameterisations, which are not necessarily physically motivated. 
These may be adequate for work with currently available data, but as redshift 
cover increases and sensitivity improves, better parameterisations will be required. 
However, it is important to highlight that, even with current sensitivity, the techniques 
described here have a notable potential for constraining specific cosmological and 
high-energy physics models beyond the standard paradigm if specific 'model-tailored' 
parameterisations are used. Within a given model, the underlying physics often points 
to specific parameterisations for the violation of standard laws, and these parameters 
can be directly related to fundamental/microphysical parameters of the underlying 
theory. This approach was initiated in reference [19], where specific parameterisations 
for distance duality violation were adopted for different models (photon-axion mixing, 
hidden photons, mini-charged particles). Other scenarios that can be probed with the 
tools described in this paper, through the use of more realistic parameterisations as 
suggested by the theory, include varying a cosmologies and dynamical dark energy. 
These will be discussed in detail in a follow-up publication [7i] . 
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